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Abstract 

Gamow solutions are used to transform self-adjoint energy operators by means 
of factorization (supersymmetric) techniques. The transformed non-hermitian op- 
erators admit a discrete real spectrum which is occasionally extended by a single 
complex eigenvalue associated to normalized eigensolutions. These new Hamilto- 
nians are not pseudo-hermitian operators and also differ from those obtained by 
means of complex-scaling transformations. As an example, Coulomb-like poten- 
tials are studied. 

1 Introduction 

Complex energies were studied for the first time in a paper of Gamow concerning the 
alpha decay (1928) [1]. In a simple picture, a given nucleus is composed in part by alpha 
particles {"^He nuclei) which interact with the rest of the nucleus via an attractive well 
(obeying the presence of nuclear forces) plus a potential barrier (due, in part, to repulsive 
electrostatic forces). The former interaction constrains the particles to be bounded while 
the second holds them inside the nucleus. The alpha particles have a small (non-zero) 
probability of tunneling to the other side of the barrier instead of remaining confined to 
the interior of the well. Outside the potential region, they have a finite lifetime. Thus, 
alpha particles in a nucleus should be represented by quasi-stationary states. For such 
states, if at time t = the probability of finding the particle inside the well is unity, in 
subsequent moments the probability will be a slowly decreasing function of time (see e.g. 
Sections 7 and 8 of reference [2]). 

In his paper of 1928, Gamow studied the escape of alpha particles from the nucleus 
via the tunnel effect. In order to describe eigenfunctions with exponentially decaying time 
evolution, Gamow introduced energy eigenfunctions i/jq belonging to complex eigenvalues 
Zg = Ec — iTG, Tg > 0. The real part of the eigenvalue was identified with the energy of 
the system and the imaginary part was associated with the inverse of the lifetime. Such 
'decaying states' were the first application of quantum theory to nuclear physics. 

Three years later, in 1931, Fock showed that the law of decay of a quasi-stationary 
state depends only on the energy distribution function uj{E) which, in turn, is meromor- 
phic [2]. According to Fock, the analytical expression of uj{E) is rather simple and has 
only two poles E = Eq ± iT, T > (see equation (8.13) of reference [2]). A close result 
was derived by Breit and Wigner in 1936. They studied the cross section of slow neutrons 



and found that the related energy distribution reaches its maximum at Er with a half- 
maximum width Fr. a resonance is supposed to take place at Er and to have "half-value 
breath" Tr [3] . The resonances can be defined as eigensolutions ■0R of the Hamiltonian 
with complex eigenvalue zr = Er — iTR/2. This complex number also corresponds to a 
first-order pole of the S matrix [4] (for more details see e.g. [5]). However, as the Hamil- 
tonian is a Hermitian operator, then (in the Hilbert space Ti.) there can be no eigenstate 
having a strict complex exponential dependence on time. In other words, decaying states 
are an approximation within the conventional quantum mechanics framework. This fact 
is usually taken to motivate the study of the rigged (equipped) Hilbert space Ti [6]. The 
mathematical structure of Ti lies on the nuclear spectral theorem introduced by Dirac in a 
heuristic form [7] and studied in formal rigor by Maurin [8] and Gelfand and Vilenkin [9]. 

Some other approaches extend the framework of quantum theory so that quasi- 
stationary states can be defined in a precise form. For example, the complex-scaling 
method [10-12] (see also [13]) embraces the transformation H — > SHS~^ — Hq, where 
S is the complex-scaling operator S — e"^''^, [r,p] = i, such that Sf{r) = /(re*^). This 
transformation converts the description of resonances by non-integrable Gamow states 
into one by square integrable states (A relevant aspect of the method is that it is pos- 
sible to construct a resolution to the identity [14]). Thus, the complex-scaled resonance 
eigenfunctions are ^-dependent so they can be normalized. Moreover, as the complex 
eigenvalues are ^-independent, the resonance phenomenon is just associated with the 
discrete part of the complex-scaled Hamiltonian [15] (but see [13]). 

In this paper we show that Gamow (decaying) eigensolutions can be used to transform 
Hermitian Hamiltonians into non-self adjoint energy operators with purely real spectrum 
or admitting a single extra complex eigenvalue with square-integrable wavefunction. The 
new Hamiltonians could be profitable as testing operators in diverse approaches includ- 
ing complex-scaling and pseudo-hermitian [16] transformations. As we shall see, it is 
not necessary to work in a equipped Hilbert space framework because the Gamow solu- 
tions will be used merely as mathematical tools. Moreover, the exponential growing of 
the Gamow solutions for large distances will be primordial in order to get well-behaved 
complex potentials. The mechanism we are going to use is the factorization method in a 
'complex' version [17]. As usual, the procedure and results can be interpreted in terms of 
supersymmetric quantum mechanics. 

The next section introduces general expressions for transforming spherically symmetric 
potentials in terms of appropriate Gamow vectors. It is shown that new complex potentials 
are derived so that their discrete spectrum is real. The Coulomb potential is managed as 
example. The Section 3 shows how the approach can be generalized to include an extra 
single complex eigenvalue into the initial discrete spectrum. The related eigensolution is 
then shown to be of finite norm. Finally, Section 4 is devoted to the concluding remarks. 
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2 Super symmetric Gamow transformations 



2.1 The complex factorization 

Let us consider the time-independent Schrodinger equation for a spherically symmetric 
potential V{r). After separation of angular variables, the equation reduces to a differential 
equation involving only the radial variable: 

H,tP{r,i)=Etlj{r,£), (1) 

which can always be integrated numerically. The reduced Hamiltonian reads 



d"^ , , d^ £(£ +1) 



where the effective potential Vi{r) has the domain Dy = [0, +oo) and the units of energy 
and coordinates have been properly chosen. 

The nature of the energy spectrum of He may be deduced from the asymptotic be- 
haviour of the solutions ipir^t) which are regular at the origin. If V{r) approaches zero 
asymptotically faster than 1/r: lim^^oo TV{r) = 0, then the energy spectrum contains two 
parts: (a) Negative discrete values Ei{i), E2{i), . . . To each of them corresponds a radial 
wavefunction of finite norm, (b) Unbound continuous positive spectrum, with solutions 
regular at the origin but indefinitely oscillating in the asymptotic region. On the other 
hand, if V{r) approaches zero as 1/r when r —>■ oo, the essential result concerning the 
nature of the spectrum persists [18]. We shall concentrate on the discrete spectrum by 
assuming that a complete set of normalized wavefunctions ipnifA) ^ has been given 
for each V{r), otherwise would not be an observable. 

We look for a complex-type factorization [17] of the Hamiltonian ([2]): 

He = AB + e (3) 

with factorization constant C 9 e = ei + i €2; ei, €2 7^ G M and a couple of not mutually 
adjoint first order operators 

A:=-^ + (3, B:=^+(3 (4) 
dr dr 

where /3 is a complex-valued function fulfilling the Riccati equation 

-P'{r) + P\r) + e = Ve{r). (5) 

This last equation is easily solved by means of the logarithmic transformation f3{r) = 
— ^ In 'u(r), with u{r) the eigensolution of He belonging to the complex eigenvalue e = 
-k^, C3 k = ki + ik2; ki,k2 G M. 

Remark that Hj = B'^A^ + e = Hi (the bar stands for complex conjugation) because 
the Hamiltonian is assumed to be self-adjoint in the Hilbert space Ti. A relevant aspect 
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of the complex factorization ([3])- ([5]) is that the reverse ordering of the factors gives rise to 
non-hermitian second order differential operators: 

BA + e = He + 2(3'{r) := he. (6) 

Conventional factorizations assume a priori A = B'^ and real e (see e.g. [19]). In coun- 
terdistinction, complex factorization is more in the spirit of the 'refined factorizations' 
reported recently [20] (see also [21]). The following intertwining relationships hold 

h^B = BHi, HeA = Ahe (7) 

which permit to determine the solutions \1/ oc B(p of he"^ = A\l/, A G C, by giving the 
solutions (fi of H£(p = \ip. The operator A reverses the action of B. In the supersym- 
metric language, and are understood as supersymmetric partners while (3{r) is the 
superpotential (see e.g. [22] and references quoted therein). 

In general, we want to keep the physical interpretation of \I' as connected with the 
probability density p(r) = |^(r) p in 7i (The dependence of \l/ on £ will be always implicitly 
considered). Hence, we look for functions 

v^oc5^=^^^ (8) 

u 

which are square-integrable in Ti (the symbol Vr(-, ■) stands for the wronskian of the 
involved functions). Of course, this last condition is not imperative in equation ([8]). For 
instance, one could extend the initial boundedness condition \ip{r,t)\'^ < oo to better 
admit another kind of normalization in order to generalize selfadjointness (e.g., in the 
picture of a equipped Hilbert space 7Y). But, in this way, the physical interpretation of 
either ?/'(r, €) or \l/(r) as wavefunctions is less clear (one dimension plane waves, for exam- 
ple, are known to be not in L^(]R) but having a probability density which is everywhere 
finite in the Dirac sense. In other words, the plane waves could be understood as energy 
Dirac vectors in Ti. However, if we apply realistic vanishing boundary conditions at a; = 
and X = L, or L-periodic boundary conditions, the plane waves can be normalized in the 
conventional form. Thus, 'free particles' are but an abstraction from the actual quantum 
world) . 

As it could be expected, the set of eigenvectors ([HD is uncommon in H: though they can 
be normalized, their elements are not mutually orthogonal [17] (An optional bi-orthogonal 
basis has been recently discussed in [23]). These vectors are natural in the spaces with 
an indefinite metric as studied in the Pontrjagin-Krein formalism [24] (see also [25]). 

2.2 Gamow transformations 

Let us show how the Gamow solutions can be used as transformation functions u{r) in 
equation ([8]). First, following Gadella-de la Madrid, we define a Gamow function as 
a solution of the time-independent Schrodinger equation with complex eigenvalue and 
purely outgoing boundary conditions [26]. Thus, if u{r) is such that u{r = 0) = 0, 
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u{r +00) ~ e-'^^ {ki < 0), and solves with E = e E C, then u{r) is a Gamow 

solution (Observe that e does not necessarily correspond to the poles of the S matrix!). In 
the context of the alpha decay, the condition u{r = 0) = describes the 'creation' of alpha 
particles inside the nucleus and obeys the fact that there cannot be any transmission into 
the region r < because the effective potential is infinite there (i.e., this condition avoids 
the incoming probabilities and is related with the adjointness of the Hamiltonian [26]). 
On the other hand, the outgoing boundary condition ensures the decay rate obeyed by 
the particles after tunneling the electrostatic barrier. 

Let us take Re(e) = Er = k'^ - kj > in e = - kf) - 2i kik2. Thus \k2\ > \ki\. We 
can distinguish two general cases: 

1) ki < 0, ^2 < 0. Here = Er — i r^/2, with = 4:kik2 > 0, is associated with the 
decaying part of the solution U{t)\(f)^-) = e"**-^^ e~*^ 

2) ki < and k2 > 0. The complex energy = Er + i T~^/2, with r+ = 4|fci| ^2 > 0, 
is associated with the growing part of the solution U{t)\(j)^+) = e~^^^^ e*"^^ . 

In both cases the roles are interchanged under complex conjugation. Now, if cor- 
respond to the poles of the S matrix, then the lifetime r = I/Fjlj decreases as the 
energy increases. Thus, for small widths (large lifetime) the energy resonances are close 
to the real axis and the Gamow vectors could be considered as bounded states for certain 
physical phenomena. On the other hand, as increases, the resonances move away from 
the real axis and the Gamow vectors are far to be considered as representative of bound 
states. 

Now, let us analyze in detail the equation ([8]). Our goal is to characterize the spectrum 
of hi as well as its eigenfunctions in terms of the analytical behaviour of f{r,i) and the 
boundary conditions of u{r). 

A direct calculation shows that u{r) oc r^+^ satisfies u{r = 0) = 0. Thereby, equation 
(ED reads 

i+1 

^(r « 1) ~ ip'(r « 1) (p{r << 1). (9) 

r 

It is clear that \I'(r) will be regular at the origin if (p is such that if{r << 1) ~ r*, s > 1. 
In other words, if (p is regular at the origin then \l/(r = 0) = 0. 

The purely outgoing boundary condition, in turn, is equivalent to the following ex- 
pression (see [26] p 630): 

lim — In ■u(r) = — lim /?(r) = —k. (10) 

r^oo (XT r^oo 

Hence, equation ([8]) reduces to 

lim \I' oc lim (f' + k lim (f. (11) 

r— >oo r—*oo r^oo 

As the solution Lp grows indefinitely as one of either e^'^'', k = \/—X, we can identify the 
following cases: 
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I) For a (denumerably infinite) set of negative discrete values A G {En{i)}, the solution 
(p in (ITT!) behaves as ~ e"""", k > 0. Thus \& ~ (A; — /t)e~'"', k > 0. 

II) If A > 0, then ~ sin(Kr — y + with 5^ the phase shift. Thus, ip is an acceptable 
eigensolution of He for any A > and represents an unbound state [18]. Hence, if 
A > then \l/(r) indefinitely oscillates when r ^ oo. 

Ill) If A G C then equation (fTTj) gives ^± ~ (±/t + k)e^'^^. Moreover, if A = e (equiv- 
alently k = k) then = and \E'+ ~ 2ke^^. The former solution is rather trivial 
as W{u,u) = in equation ([8]). Now, as ki < 0, it seems that could satisfy 
limr-->+oo l^+l = 0. However, in such a case, (p+ should also satisfy both conditions 
<^+(0) = and ip+ oc e''^, ki < 0, which is not possible since A is complex and He 
is a selfadjoint operator in Hep>+ = A(y9+. A similar situation arises for any complex 
number A different from e. 

In summary, for Gamow transformation functions in ([H]), if G L^(M"'") then \I' G 
L^(M+). Furthermore, equation (jSj) does not produce eigenfunctions of the non-hermitian 
Hamiltonian he belonging to complex eigenvalues. Thereby, the complete discrete spec- 
trum ad{He) of the initial Hamiltonian He is inherited to the Gamow transformed Hamil- 
tonian he. In order to exhaust our analysis, let us consider the complex factorization ([6]). 
It is easy to verify that the kernel of A provides an eigenfunction C,t{r) of he belonging to 
e G C. Thus, oc l/u fulfills he^^ = ^(.t- However, as u is a Gamow vector, diverges 
at the origin as r^^"^. In other words, is out of H and it is deprived of a physical 
meaning. The same situation arises by considering the two-dimensional kernel of the 
product BA. Hence, there are no more square-integrable solutions of he and the discrete 
spectrum adihe) is just the same as ad{He). 

We have then constructed a non-hermitian Hamiltonian he which is strictly isospectral 
to the initial spherically symmetric Hermitian Hamiltonian He- A simple calculation shows 
that the global behaviour of the new potential veir) = Ve{r) — 2u'{r)/u{r) is as follows 



Thus, for small distances, a particle with energy En{i) interacts with the field as having a 
quantum number In the asymptotic region the particle behaves as free of interaction. 
On the other hand, the intermediate region could be interpreted as 'opaque' in the sense 
that the particle interacts with a series of wells and barriers which alternate their positions 
in the real and imaginary parts of ve{r) (see the discussion on the optical bench given 
in [27]). The next section elucidates the applications of the method by transforming the 
Coulomb potential. 

2.3 Non-hermitian Hamiltonians with hydrogen— Uke spectrum 

If the radial potential in ([2]) is the Coulomb one V{r) = — 2/r, the convenient Gamow 
vectors are given by the expression (see Figured]): 




(12) 



= /+^e-^"iFi(£ + 1 - 1/fc, 2£ + 2, 2kr) 



(13) 
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with i-Fi(a, c, z) the Kummer's function. The units of energy and coordinates are respec- 
tively taken as £ = Ze^/2rB (= Z 13.5 ev) and = h^/Ze^m (= 0.529 x IQ'^/Z cm). 
The solutions (fT3|) have been explicitly derived in [17]. 




-2-10 1 2 3 



Figure 1: The Argand-Wessel diagram of the hydrogen's Gamow vector u(r), from zero to 20 Bohr radii 
tb, with £ = 1 and e = —0.2604 + i 0.104 (x 13.5ev). Horizontal scale stands for the real part. The disk 
is at r = 1 and the circle dX r — 19 r^. 

Once these Gamow vectors have been used as transformation functions in ([H]), the 
non-hermitian potential v^{r) resembles a cardiod curve as depicted in the complex plane 
(see Figure [2]). Notice that v^{r) becomes almost real for small distances and goes to 
+00 on the real branch for r = 0. On the other hand, this potential goes to zero when 
r +00. The imaginary part of vi{r) becomes relevant for intermediate distances (i.e. 
at distances which are between 2 and 6 Bohr radii, for the parameters considered in the 
figure) . 
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Figure 2: The Argand-Wessel diagram of the non-hermitian potential for the same parameters 
as in Figure [1] This potential has a discrete spectrum identical with that of the hydrogen atom. The 
right-hand side figure is a detail of the cardiod-type one. The disk is at r = 2r_B and the circle at 
r = 6 tb. 

Finally, Figure [3] depicts the potential Vf+i(r) as well as the real part of f^(r). Observe 
the presence of barriers and wells in the intermediate distances. A similar situation occurs 
for the imaginary part of Vi{r). These 'partial potentials' induce local 'resonance' effects 
which are not present in the Hermitian potential V^+i(r). Thus, the spatial distribution 
of the wave-packets corresponding to Vi{r) differ from that of the wave-packets of V£+i(r) 
at the same energy. 
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Figure 3: The real part (continuous curve) of the cardiod-type potential in Figure [2] contrasted with 
the effective Coulomb potential Vi=2(r). 

3 Generalized Gamow transformations 

As it has been shown in the precedent section, though the Gamow vectors u{r) could 
have a definite physical meaning as resonant states of Hg, we consider them merely as 
transformation functions to construct the non-hermitian Hamiltonians hg. In general, all 
the unphysical (not square-integrable) solutions of the Schrodinger equation are useful to 
construct new Hamiltonians admitting real spectra and square-integrable eigenfunctions 
[28,29]. In particular, if the factorization constant e is a real number, the conventional 
factorization operators A = B"^ are automatically recovered. 

Now, we extend the previous results by opening the chance to incorporate complex 
eigenvalues with square-integrable solutions in the spectra of the transformed Hamil- 
tonians. First, notice that the general solution of ([I]) is, for small distances, a linear 
combination of two particular solutions: r^~^^ and r~^. The second one is usually rejected 
because it is singular. Moreover, in the context of alpha decay, a vector u{r = 0) 7^ does 
not describe the 'creation' of alpha particles. We shall relax the Gamow condition at the 
origin to include the solution but preserving the purely outgoing condition e~^^. Let 
us remark that a 'generalized' Gamow vector Lv{r), satisfying these new conditions, still 
is unphysical in the sense that it is not square-integrable in Ti. 

The relevance of our generalization lies on the fact that expressions ([I])-(IH]) still hold if 
uj{r) is taken as the transformation function. Equation ([9]), on the other hand, is slightly 
modified: 



Hence, the same conclusion is obtained: if (/? G L^(M) then \Ef G L^(]R). However, for com- 
plex eigenvalues of h£{r), the kernel of A provides the eigensolution oc l/o;, which can 
be normalized in Ti. and satisfies hi'^^ = e \E'e. It is easy to check that the corresponding 
complex conjugate is neither in the kernel of A nor that of BA (see equation ([6])). In 
counterdistinction, if ip{r) is eigensolution of Hg belonging to e, then ip{r) belongs to e as 
Hi is selfadjoint. 

Therefore, the discrete spectrum of hi is now given by ad{Hi) U {e}. On the other 
hand, the new potential Vi{r) = Ve{r) — 2uj'{r)/uj{r) behaves in this case as 




(14) 



r 




(15) 
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with a similar interpretation as for (fT2|) but changing i + 1 hj £ — 1. 




(a) (b) 

Figure 4: (a) The generalized hydrogen's Gamow vector uj{r), plotted from zero to 20 Bohr radii and 
the same parameters as in Figure [T] The disk is at r = 0.05 and the circle at r = 19.5 r^. (b) The 
square-integrable wavefunction belonging to the complex eigenvalue e = —0.2604 + i 0.104 (xl3.5ev). 
The disk is at r = 0.05 and the circle at r = 19 rg. 

Figures H] and [5] show respectively the behaviour of the generalized Gamow vector 
ui{r), the wavefunction "^eif) and the new non-hermitian potential f^(r) for the Coulomb 
case V{r) = —2/r. The related transformation function is [17]: 

cu(r) = r^+^e-''''[iFi{i+l-l/k,2i + 2,2kr) 

(16) 

+ ^U{i + l-l/k,2i + 2,2kr)] 
where ^ is a complex constant and U{a, c, z) is the logarithmic hypergeometric function. 
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(a) (b) 

Figure 5: (a) The non-hermitian potential constructed via the generalized Gamow vector of 
Figure S^. Here we have plotted from zero to 40 Bohr radii. The disk is at r = 2rs and the circle at 
r — Grs- This potential has the same spectrum as the hydrogen atom extended by an extra complex 
eigenvalue at e with the square-integrable wavefunction of Figure |4)3. (b) The corresponding real part 
(continuous curve) contrasted with the Coulomb potential Vf=o('')- 



4 Concluding remarks 

The Gamow (decaying) eigensolutions have been shown to be appropriate transforma- 
tion functions in the framework of supersymmetric quantum mechanics. Non-hermitian 
Hamiltonians, which are supersymmetric partners of spherically symmetric self-adjoint 
energy operators, have been constructed so that they admit purely real spectrum with nor- 
malized wavefunctions. For other of these new non-hermitian operators an extra complex 
eigenvalue with square-integrable eigensolution is present. 
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At first sight, our results could be connected with those derived in the complex-scaling 
method. However, the Gamow-transformed potential ve{r) = Vi{r) — 2uj'{r)/uj{r) is not 
as simple as the complex-scaled potential Ve — e*^^ V^(re'^). In general, an intertwined 
Hamiltonian hB = BH, HA = Ah (factorized in a refined way [20]: H = AB + e, 
h = BA + e) could correspond to a complex-scaled Hamiltonian hgS = SH if hg = h. 
Thus, there must exist a couple of differential operators M = AS and = BS~^, such 
that [H, M] = [h, iV] = 0. A particular case has been recently reported [30] (see also [31]) 
by considering conventional factorization operators = "F^ + Q;(r), a : R i— > IR, real 
factorization constants and the squeezing operator S — U^ — e'5{^'P}^ [r,p] = i. The so 
derived 'scaled intertwined' Hamiltonian hx has a real potential v\{r) = e^^ V{e^r) —a'{r) 
and real discrete spectrum (Td{h\) = {e^^E, e^^i?„}„gN, where e^^£ is the ground state 
energy and G aa{H). This procedure allows to deform the excited energy levels of 
hx but leaving unaffected the ground state Eq: ad{hx) ^ {£q, {£o/£)En}. This is a 
remarkable profile of the factorization which is rarely considered in the literature. Thus, 
it seems that complex-scaled Hamiltonians he could be successfully constructed as an 
application of the technique reported in [30]. Work in this direction is in progress. 

On the other hand, the possible connection of our results with other approaches as 
the PT-symmetry [32] or the pseudo-hermitian transformation [16] has been discussed 
in a previous work [17]. 

Finally, we have presented the case of first order, not mutually adjoint, intertwining 
operators A, B. However, the method can be iterated at will by considering h^ as the new 
initial Hamiltonian. The nth iterated result can be also obtained by means of intertwining 
operators of nth order. In particular, the second order case can be properly used to obtain 
self-adjoint Hamiltonians with spectrum identical to the initial one [17,33]. It is also 
possible to show that second order transformations can produce non-hermitian operators 
with real spectrum extended by two extra complex eigenvalues e and e. These results will 
be published elsewhere. 
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